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ABSTRACT 

The densest infrared dark clouds (IRDCs) may represent the earliest observable stage of high-mass star formation. These clouds are 
very cold, hence they emit mainly at far-infrared and sub-mm wavelengths. For the first time, Herschel has provided multi- wavelength, 
spatially resolved observations of cores within IRDCs, which, when combined with radiative transfer modelling, can constrain their 
properties, such as mass, density profile and dust temperature. We use a 3D, multi-wavelength Monte Carlo radiative transfer code 
to model in detail the emission from six cores in three typical IRDCs seen in the Hi-GAL survey (G030.50+00.95, G03 1.03+00.26 
and G03 1.03 +00.76), and thereby to determine the properties of these cores and compare them with their low-mass equivalents. We 
found masses ranging from 90 to 290 M© with temperatures from 8 to 11 K at the centre of each core and 18 to 28 K at the surface. 
The maximum luminosity of an embedded star within each core was calculated, and we rule out the possibility of significant high 
mass star formation having yet occurred in three of our cores. 

Key words, stars: formation - ISM: clouds - ISM: dust 



1. Introduction 

We present observations, performed with the ESA Herschel 
Space Observatory ([Pilbratt et al.||201Q|), of a section of the 



Galactic Plane. We use HerscheVs large collecting area and 
powerful science payload to perform imaging photometry us- 
ing the PACS (P hotodetector and Array Camera Spectrometer; 
[Poglitsch et~ari[2010^) and SPIRE (S pectral and Photometric 
Imaging Receiver; Griffin et al. 2010| ) instruments. These ob- 
servations were carried out as part of HerscheFs Science 
Demonstration Phase (SDP) for the Herschel Infrared Galactic 
Plane Survey (Hi-GAL; Molinari et al. 2010a b ) open time Key 
Project to map the inner Galactic Plane. In particular, we are 
studying cores within the infrared dark clouds (IRDCs) found in 

the Galactic Plane. 

IR DCs were intially disco vered by the ISO ([Perault et al. 
T996] ) and MSX ( [Carey et al. 1998; E gan et al.||1998| ) surveys 
as regions of high contrast against the mid-infrared (MIR) back- 
ground. The densest IRDCs may eventually form massive stars, 
and are therefore presumed to represent the earliest observable 
stage of high mass star formation. Some IRDCs contain cold, 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



compact cores, which are believed to be the high mass equiva- 
lent of low-mass prestellar cores. 

IRDCs and their cores have low temperatures (T<25 K; e.g. 
Egan et al.|1998[ Teyssier et al.|2QQ2 ) and therefore emit mostly 
at far-infrared (FIR) and sub-nmi wavelengths. However, until 
now, spatially resolved observations have not been available in 
the FIR. Herschel changes this, and offers high resolution FIR 
images of not only IRDCs, but of the individual cores within - 
potentially allowing, for the first time, temperatures to be deter- 
mined and column densities and masses to be accurately con- 
strained for a large number of IRDCs. In this paper we present 
some first examples of how the Herschel data can be used to 
trace the physical parameters of cores within IRDCs. 



2. A Sample of IRDC Cores in the Hi-GAL Science 
Demonstration Fields 

Hi-GAL is an Open Time Key Project of the Herschel Space 
Observatory and is performing a survey of the inner Galactic 
Plane (|/|<60°, \b\<r) using the PACS and SPIRE instruments. 
The two are used in parallel mode to map the Milky Way Galaxy 
simultaneously at five wavelengths, from 70 to 500 yum, with res- 
olutions up to 5'' at 70 yum. 
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Table 1. The physical properties of the six cores. 







G030.50+00.95 


G03 1.03+00.26 


G03 1.03+00.76 






Core A 


CoreB 


Core A 


CoreB 


Core A 


CoreB 


R.A. (2000) 




18:43:34.7 


18:43:34.9 


18:47:01.5 


18:47:04.2 


18:45:10.2 


18:45:17.3 


Dec. (2000) 




-01:42:29 


-01:44:36 


-01:34:36 


-01:33:39 


-01:20:00 


-01:19:12 


Name in PF09 




SDC30.442+0.958 


SDC3 1.039+0.241 


Not included 




Distance (pc) 




3400 


3400 


4900 


4900 


3400 


3400 


Semi-major Axis (pc) 




0.6 


0.7 


0.45 


0.55 


0.4 


1.0 


Ro (pc) 




U.Uo 


U.U/ 


U.U4 


U.Uj 


0.04 


0.1 


Aspect ratio 




1.93 


3.16 


2.46 


2.22 


1.70 


1.65 


Position Angle 




50 


85 


125 


340 


335 


85 


^o(H2) (cm-3) 




1.4x10^ 


4.5x10^ 


3.8x10^ 


4.2x10^ 


6.0x10^ 


2.5x10"^ 


ISRFnabing 




8.5 


18 


18 


81 


1 n 
1 / 


A Q 
4.0 


^6=90°/'^ 9=0° 




1.93 


3.16 


2.46 


2.22 


1.59 


1.65 


^ 500 lum K-fy) 




5.0+0.5 


9.0+0.9 


3.3+0.3 


10+1 


4.5+0.5 


4.4+0.4 


Model Temperature (K) 




9-19 


11-22 


9-22 


11-28 


8-22 


10-18 


SED Temperature (K) 




14+1 


16+2 


14+1 


17+2 


14+1 


14+1 


Model mass (M©) 




120 


140 


170 


290 


110 


90 


SED mass - /c=0.05 cm^ ^ 


f' (Mo) 


130 


170 


170 


310 


110 


100 


SED mass - /c=0.03 cm^ | 


f' (Mo) 


210 


280 


280 


510 


180 


160 






0.75 




1.50 




1.25 






<6 


<3 


<200 


<1000 


<40 


<9 



Notes. Rq is the flattening radius within the core, ^o(H2) is the central d ensity of each core, l^RFuabing is the ISRF of the core at FUV given in 
multiples of the Habing flux — i.e. 1.6x10"^ erg s~^ cm~^ ( [Habing 1968| ), Tq^^^o /t^^qO is the ratio of optical depths at ^=90° and at ^=0°and Fsooixm 
is the flux density of each core at 500 iim. The model temperature is the range of temperatures in the core calculated from the model, a nd the SED 
temperature was calculated from single-temperature greybody fitting. The model mass was calculated using /c=0.05 cm^g~^ at 500 yum (Ossenkopf 
[& Henning|1994| . The two SED masses were calculated using the method of ,Hildebrand| |1983), th e SE D temperature and the flux density at 



500 jjm (see Section 
star that could be emi 



The diff'erence in ISRFs, ln(j-h calculated from the model (see Section 
"ded in each core (see Section l43l. 



4.3 



Ljnax is the maximum luminosity of a 



Two regions of the Hi-GAL survey were completed during 
SDP time. The regions are each four degrees square and center 
around /=3Q°, /7=0°and /=5 9°, b=Q\ Using th ese observations 
and th e |Simon et al.| ( |2Q06al hereafter 806) and |Peretto & Fuller] 
(2009 hereafter PF09) catalogues, IRDCs can be identified in 
the Galactic Plane with none of the uncertainty that exists in 
using MIR observations alone (e.g. being unable to diflPerentiate 
between an IRDC and a dip in the MIR background). 

PACS data reduction was performed using the Herschel 
Interactive Pipeline Environment (HIPE; Ott 20 1()|), albeit with 



some departures from the standard processing ( Poglitsch et al. 
2010| ). Most notably, the standard deglitching and crosstalk cor- 



rection were not used due to poor results in these fields, and 
custom procedures were written for drift removal. 

SPIRE data processing required less devi ation from the stan- 
dard processing methods ( GriflSn et al.||201 0), with both stan- 
dard deglitching and drift removal producing acceptable results. 
In both case s, the RO MAGAL Generalised Least Squares algo- 
rithm (T raficante et al.^2010) was used to produce the final maps. 
A more thorough discussion of the entire data reduction process 
can be found in " Traficante et al.| ( |2010| ). 

Each of the S06 sources in these regions was studied at 6 
wavebands, i.e. S pitzer GLIMPSE Syum ( Benjamin et al.||2003[ 
Werner et al.|2004 ) and the five Hi-GAL wavelengths (70, 160, 
250, 350 and 500 yum). If the source showed extinction at MIR 
but emission in the FIR it was identified as an IRDC. Within 
the /=30° region, approximately 330 (PF09) IRDCs were found, 
each containing at least one infrared dark core. 

Our aim is to model those cores without embed- 
ded protostars. We focus on three objects, G030. 50+00.95, 
G03 1.03+00.26 and G03 1.03+00.76, which meet the criteria for 
IRDCs. All three contain two cores each, designated 'Core A' 



and 'Core B' in order of Right Ascension within each cloud. 
Positions and physical properties of all six cores can be found in 
Table [T] and images are displayed in Figs.[T}|3] 



2.1. G030.50+00.95 

In this cloud. Core B is the more extended and brighter of the 
two cores and located to the south of Core A (see Fig.[T]). In the 
250yum and 350yum wavebands, it is possible to see that Core 
B can be separated into two cores. However, as the individual 
peaks cannot be seen at either 160yum or 500 yum, the two are 
grouped together and modelled as a sing le core. A distance o f 
3.4 kpc for the entire cloud is assumed ( Simon et al.||2QQ6b 



No previously calculated masses are available for these cores. 
Results calculated from our data are presented in Table [T] 

2.2. G031 .03+00.26 

S06 find three cores in this cloud. However, only two can be seen 
as significant emission sources in the FIR and so here the third 
is ignored. Core B is the more extended and brighter of the two 
and is located to the north of Core A (see Fig. [2]). Emission from 
Core B can be clearly seen at wavelengths as short as 160yum, 
whereas emission from Core A is not seen shortward of 250 jim. 



_We take a distance of 4.9 kpc for both cores ( [Teyssier et al. 



2002j) although we note that there are two emission line com- 
ponents among the line of sight, with the second emission line 
giving a distance of 6.6 kpc. There fore there is some uncer- 
tainty in this distance. Parsons et al.| (|2QQ9 ) use 850yum data and 
find masses of 420 and 560 M©, respectively, for these cores. 
However, these authors used a diflTerent size of aperture. We 
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L. A. Wilcock et al.: A Radiative Transfer Model of IRDCs in the Hi-GAL Survey 
Table 2. The observed, integrated flux densities of the six cores. 



Wavelength 






Flux Density (Jy) 








G030.50+00.95 


G03 1.03+00.26 


G03 1.03+00.76 




Core A 


CoreB 


Core A 


CoreB 


Core A 


CoreB 


70yL/m 


<0.03 


<3.6 


<0.5 


<10.1 


<1.0 


<1.3 


160yt/m 


13.4 


41.9 


10.9 


60.3 


13.7 


15.4 


250yL/m 


18.0 


31.4 


11.7 


47.5 


13.4 


13.3 


350yt/m 


9.6 


16.0 


6.8 


23.0 


8.8 


8.8 


500yL/m 


5.0 


9.0 


3.2 


9.7 


4.5 


4.4 



make a direct comparison with these results in Section 4.1 be- 
low. Results calculated from our data are presented in Table T] 



2.3. G031 .03+00.76 

Here, Core B is located to the north of Core A (Fig. [Sj and, 
although it is the more extended of the two, it is actually fainter 
than Core A. Neither core can be seen in emission at wavelengths 
shorter than 160yum. A distance of 3.4 kpc is assumed (Simon^ 
et al.|2006b| ). No previously calculated masses are available for 



these cores. 

A third core in this cloud, located to the south of Core A, is 
not modelled here due to a bright Syum source at the same posi- 
tion, implying the core may have an internal heating source. A 
possible fourth core exists to the north of Core B, we ignore this 
as S06 state that it belongs to a difl'erent IRDC with no known 
distance. Results calculated from our data are presented in Table 

m 



3. IRDC Models 

We model the cores in two difl'erent ways: first, using a greybody 
single-temperature fit, and second, using the 3D Monte Carlo ra- 
diative transfer code, Phaethon ( [Stamatellos & Whitworth|2003[ 
[20051 IStamatellos et al.,201Q^ 

3.1. Single Temperature Fitting of SEDs 

The flux density, integrated over twice the FWHM of each core 
using an elliptical aperture, was measured in all five FIR maps 
(see Table [2]) and an SED (spectral energy distribution) was 
plotted. These flux densities have been background subtracted, 
where the background was defined using an off'-cloud, elliptical 
aperture. For example, with G03 1.03+00.76 the area to the right 
of Core B was used. The background subtraction removed up 
to 50% of the original flux, the exception being at 70 yum where 
over 90% of the original flux was removed. A single-temperature 
greybody (modelli ng thermal emissi on from cold dust) was fit- 
ted (using MPFit; IMarkwardt||2008 ) to each core, shown as a 
in Fig. [4j ■ ' ' 



dashed line i 



This has the form 



(1) 



where: Fy is the flux density at frequency y; By(T) is the black- 
body function at temperature T; Q is the solid angle subtended 
at the observer by the source; Imn is the mass of a hydrogen 
molecule; X is the mass fraction of hydrogen and Ky is the dust 
mass opacity (e.g. |Kirk et al.|2QlQ[ [Ward-Thompson et al.|2010| . 
Ky is given by 



where jS, the dust emissivity index, was set to 1.85 ( |Ossenkopf| 
&Henning|r994l ). 

The model was fitted between 160yum and 500 yum, with the 
flux density at 70 jim being used as an upper limit, as the cores 
are not visible at this wavelength. The temperature was allowed 
to vary over a range 3-40 K. The best fit temperatures for the 
cores are given in Table [T] 

3.2. Radiative Transfer Modelling of the Cores 

The cores were modelled using Phaethon, a 3D Monte Carlo 
radiative transfer code. The code uses luminosity packets to rep- 
resent the ambient radiation field in the system. These packets 
are injected into the system where they interact (are absorbed, 
re-emitted or scattered) with it stochastically. The ambient ra- 
diation field is taken to be a multiple of a modified version of 
the|Black|( |1994[ ) interstellar radiation field, which gives a good 
approximation to the radiation field in the solar neighbourhood. 

The input variables of the code are the strength of the ambi- 
ent radiation field, the density profile, the size and geometry of 
the core (i.e. spherical, flattened or cometary) and the dust prop- 
erties of the system. The code calculates the temperature profile 
of the system as well as SEDs and intensity maps, at diff'erent 
wavelengths and viewing angles. 

All six cores showed some measure of eccentricity in the 
observations and so we re modelled with a flattened geometry — 
see Stamatellos et al.| ( ^2004^ ,2Q1QJ for details. In this case the 
density profile is given by: 



n(r,e) = no (H2) 



1+A(^f [sin(^)F 



(3) 



Ky = OA cm-^g-^ X (y/lOOOGHz)^, 



(2) 



where r is the radial distance, 6 is the polar angle and Rq is the 
flattening radius (i.e. the radial distance for which the central 
density is approximately constant). ^o(H2) is the central density, 
which is controlled as an input variable. A is a factor that controls 
the equatorial to polar optical depth ratio and determines how 
flattened the core is. p determines how quickly the optical depth 
changes from equator to pole, and was set to 2. 

The FWHM of the major axis of each observed core was 
measured at 250 yum and used as the model core's semi-major 
axis. The flattening radius, Rq, was set at one tenth of this 
value. The central density, ^o(H2), and interstellar radiation field 
(ISRF) incident on the core were then varied until the output 
model's SED matched the observed data. The model SED is 
shown as a dotted line in Fig.|4] The final values for ^o(H2) and 
the incident ISRF are noted in Table [T] 

It should be noted that the internal structure of cores within 
infrared dark clouds have not yet been observed in detail. 
Hence our assumption of an elliptical geometry may be an over- 
simplification. If there is structure on smaller scales than we can 
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Fig. 4. SEDs of all six cores. Upper row (left-right): G030.50+00.95, G03 1 .03+00.26 and G03 1 .03+00.76 Core A. Lower row (left- 
right): G030.50+00.95, G031.03+00.26 and G031.03+00.76 Core B. The dotted line shows the SED calculated from Phaethon, 
using yS=1.85. The dashed line shows the single-temperature grey-body, with temperatures as noted in Table[T] also using yS=1.85. 



resolve this might affect our results. For example, small scale 
fragmentation might allow the ISRF to penetrate further into the 
cores, meaning that our calculated values are probably upper val- 
ues. 



Figs.[T}|3 show the output of the model at wavelengths corre- 
sponding to the wavelengths of the observed data. The modelled 
images have pixels of 0.02x0.02 pc in size, which corresponds to 
8'' pixels for G03 1.03+00.26 and W pixels for G030.50+00.95 
and G03 1.03+00.76. All the images have been convolved with 
the telescope beam. For wavelengths where no emission is vis- 
ible in the model, an image showing background radiation is 
shown. As no attempt was made to correctly model the surround- 
ing area, these images are not an accurate representation of the 
parent cloud but do show the cores in absorption as seen in the 
observations. 



With the possible exception of G3 1 .03+00.26 Core B (which 
may have evidence of an 8yum point source), the 8yum im- 
ages shows no emission in either the model or the observed 
data, as IRDCs do not emit significantly in the MIR due to 
their low temperatures. The modelled images show emission 
starting at 250yum for G031.03+00.26 Core A and 160 yum for 
the remaining five cores (albeit only faintly for the cores in 
G03 1.03+00.76) and absorption at shorter wavelengths. The 
modelled emission fits well with the observations, at all wave- 
lengths. 



4. Discussion 

4.1. Core Masses 

The masses of the cores were calculated in two ways: (i) using 
the radiative transfer models; and (ii) using the observed flux 
densities at 500 yum. The latter was carried out using the equation 



M = 



FyD^ 

By(T)Ky' 



(4) 



(Hildebrand 1983 ), where D is the source distance. We call the 
former the model mass and the latter the SED mass. 

Two difl'erent val ues of opacity at 500 //m were 
used: a:=0.03 c ni^ ([Os senkopf & Henning 1994 ) and 
K=0.05cm^g-^ ([Preibisctr et al. 1993; Henning et al. 1995). 



The latter is the value used to calculate the masses of low-mass 
prestellar cores (e.g. 'Kirk et al. '2005V The model masses were 
found using the Ossenkopf & Henning (1994 ) opacities. 

The model masses match well with the SED masses calcu- 
lated using a:=0.05 cm^ g"^ but not with those calculated using 
/c=0.03 cm^ g"^ . Both sets of SED masses and the model masses 
are shown in Table [T] 

The masses calculated here are in agreement with those 



20011 



ar, m 



found in previous studies of IRDCs (e. g. Williams et al.| 
Simon et al. 2006b{ Rathborne et al. 12006*). In particu 
Rathborne et al.| ( |2006| ), masses of dark cores were calculated 
and found to vary from 10 to 2IOOM0. All of our cores fit within 
this range. 

850 yu m masses for both core s in G03 1.03+00.26 were cal- 
culated by Parsons et al. ( |2009 ). Core A was found to have a 
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0.01 0.1 1 0.01 0.1 1 

Radius (pc) Radius (pc) 

Fig. 5. Dust temperature profiles of the cores at two diff'erent directions of the cloud, as calculated from the radiative transfer model. 
The lower line corresponds to the midplane of the flattened structure (0 = 90°), whereas the upper line corresponds to the direction 
perpendicular to the midplane (0 = 0°). Upper row: G030.50+00.95 Cores A and B; Middle row: G3 1.03+00.26 Cores A and B; 
Lower row: G3 1.03+00.76 Cores A and B 



mass of 460 Mn and Core B of 520 M^. W hen using the same 

mm and K=0.05 cm^ g-\ 



Parsons et al. 



core boundaries as in 
we find 500 yum masses of 470 M© and 540 M©, similar to the 
850//m masses. This gives us additional confidence in our cal- 
culations. 



4.2. Core Temperatures 

The energy deposited in the gas due to cosmic ray ionisation 
heating is much lower than the energ y absorbed by th e dust due 
to external heating from the ISRF (Evans et a l.|2001 ). Thus, we 
can safely assume that the dust temperature is not aff'ected by 
the energy exchange between gas and dust. The lack of MIR 
point sources also implies that heating via compression is mini- 
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mal ( B attersby et al . |20 10| ) . These cores (again with the possible 
exception of G03 1.03+00.26 Core B) therefore have no signifi- 
cant internal heating sources and are heated mainly by the exter- 
nal radiation field. 

The temperature profiles of the modelled cores are shown in 
Fig [5] They show a low temperature throughout the centre of 
the core, before it eventually begins rising, peaking at the outer 
edges. A higher temperature at the edges than at the centre is 
expected due to the lack of any internal heating source. Some 
stepping can be seen in the graphs, this is due to the temperature 
table having a spacing of 0.2 K. 

This gradient is cons istent with previous studies of IRDCs 
(e.g. Peretto et al.||2010| ) and is similar to w hat is seen in low- 
mass prestellar cores (e.g. [Nutter et al.|2008| ). However, IRDCs 
have higher outer temperatures than low-mass prestellar cores, 
which can be explained by a higher ISRF around these IRDCs 
relative to nearby regions of star formation. 

As all six cores have a flattened geometry, the temperature- 
radius profile varies with the direction inside the core. The two 
extreme cases, ^ = 0° and 6 = 90°, are shown in the temperature- 
radius profiles. 6 = 90° is along the midplane of the structure, 
while ^ = 0° corresponds to the short axis of a flattened core, 
perpendicular to the midplane. 

The average temperature from the model of each core 
agrees with the result from the single-temperature greybody fit. 
However, the temperature-radius profile shows this fit to be an 
over-simplification, with temperatures within each core varying 
by ~10K and being dependent upon angular direction in the 



core. 



4.3. The Interstellar Radiation Field 



Five of the six cores required a higher ISRF than the solar neigh- 
bourhood value in order to correctly match the observed flux 
densities at shorter wavelengths (i.e. 160 and 250yum). As an 
example we show the cores' ISRF at far-ultraviolet (FUV) in 
Table [T] This is merely to indicate how the ISRF varies from 
core to core. For comparison, th e solar neighb ourhood ISRF at 
FUV is 7.6x10"^ erg s"^ cm'^ /Black" ^199?) or 4.6 times the 
Habing flux (1.6x10"^ erg s"^ cm'^ Habing 1968 ). The IRDCs 
typically have an ISRF around four times greater, consistent with 
them being in the stronger interstellar radiation field of the inner 
Galactic Plane. 

Diff'erences between the ISRF required for cores within the 
same field can be explained by one core being more deeply 
buried within the parent cloud and thus having less external ra- 
diation heating it. The diff'erence between the ISRFs of the cores 
are listed in Table jlj as (^), where h corresponds to the core 
that is modelled using a higher ISRF and I2 to the core modelled 
with the lower ISRF (within the same field). This ratio can be 
thought as a measure of how more deeply a core is embedded 
when compared with the other core in the same field. 

Core B in G03 1.03+00.26 appears to reside in an unusually 
intense radiation field, 17 times that of the solar neighbourhood 
and far higher than any of the other cores, possibly pointing to an 
abnormally high ISRF in this field. Core A in the same cloud has 
an ISRF value comparable to the other cores, which implies that 
it is either more deeply buried within the IRDC than Core B or 
there is a highly asymmetrical external radiation field incident on 
this cloud. Alternatively, diff'erences in emission from the cores 
could be the result of diff'erent physical properties (e.g. dust grain 
size) between diff'erent cores in the same IRDC. Core B may not, 
in fact, be completely externally heated. This one is unusual in 



our sample, in that it may have an embedded protostar. There is 
some evidence for an 8 jim source associated with this core. 

4.4. Internal Heating Sources 

A lack of significant internal heating in the modelled cores was 
assumed due to the absence of any near-infrared or MIR emis- 
sion. This is further supported by the temperatures of the cores 
decreasing towards their centres. However, more detailed study 
into the heating of infrared dark cores needs to be performed be- 
fore this can be conclusively proven. For now we can put limits 
on the luminosity of any embedded source. 

Using the column-density estimated for each core, we can 
calculate the visual extinction, the K-band extinction, and finally 
the extinction at 8 yum (Ag; Inde betouw et al.|2005| ). A peak 8 jim 
flux was measured for each core (Fg) and, using 



F = FslO 



A8/2.5 



(5) 



the maximum flux density of any embedded source (F) was cal- 
culated. This was then converted to a luminosity (L) by multiply- 
ing by the Spitzer 8 fim bandwidth to convert from flux density 
to flux and using 



L = 4nD^F. 



(6) 



These values are listed in Table \T\ We note that 
G3 1.03+00.76 Core A has some 8yum sources around it's edge, 
but they do not appear to be centred on the core. 

Our upper limits for G30.50+00.95 Cores A and B and 
G3 1.03+00.76 Core B are fairly low. We can thus rule out 
any significant massive star formation having occurred in their 
centres, making it possible that these cores are the high-mass 
equivalents of low-mass prestellar cores. For both cores in 
G03 1.03 +00.26 and in G3 1.03 +00.76 Core A, the upper limits 
are too high to reach this conclusion. 

In the case of G03 1.03+00.26 Core B a very faint source can 
be seen at 8 jim from which F^ was measured. These calcula- 
tions make the assumption that this source is deeply embedded 
within the cloud. Should the source instead be in the foreground 
or towards the nearer edge of the IRDC then the maximum em- 
bedded luminosity would be very much lower. 

5. Summary and Conclusions 

We have confirmed the status of G030.50+00.95, 
G03 1.03+00.26 and G03 1.03+00.76 as IRDCs using mid- 
and far-infrared observations. Six infrared dark cores have been 
found in the FIR in these IRDCs and modelled using the 3D 
Monte Carlo radiative transfer code Phaethon. The radiative 
transfer models reproduce the observed SEDs and the FIR 
images of these cores. 

The Phaethon code calculated masses for the cores ranging 
from 90 to 290 M©. The masses calculated from the SED tem- 
perature ranged from 100 M© to 3IOM0 (with a:=0.05 cm^ g"^). 
These are in agreement with previous masses found for infrared 
dark cores. It should be noted that objects of this size and mass 
are not likely to create a single star, but rather a cluster of stars. 

The greybody fitted SED temperatures varied from 14 to 
17 K. The average temperatures calculated from the model 
agreed with these values. However, the model showed that the 
single-temperature fit was an over- simplification of the temper- 
atures inside the cores as temperature-radius profiles showed 
changes in temperature of ~ 10 to 17 K within a single core. The 
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model calculated temperatures of 8 to 1 1 K at the centre of each 
core and 18 to 28 K at the surface. 

Significant temperature gradients were also seen by |Peretto 
et al. pQ10| ), where temperature maps of 22 IRDCs were cre- 
ated and non-uniform temperatures ranging from 10 to 22 K 
were found. [Stamatellos et al.| ( |2010| ) used Phaethon to model 
an IRDC core and found temperatures from 10 to 21 K, consis- 
tent with our results. 

The maximum luminosity of a star embedded within each 
core was found. The upper luminosity limits for three of the 
cores (G30.50-hOO.95 Cores A and B and G3 1.03-^00.76 Core 
B) are low, ruling out the possibility of any significant high mass 
star formation having yet occured and making it likely that these 
cores are the high-mass equivalent of prestellar cores. 

In most cases, the amount of radiation incident on the cores 
was found to be higher than that in the local neighbourhood, 
as expected for objects within active star-forming regions in the 
Galactic Plane. Diff'erences in the amount of radiation falling on 
cores in the same cloud were attributed to one core being more 
deeply buried within the IRDC than the other. This method can 
now be used on all cores within IRDCs found in the Hi-GAL 
survey. 
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